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ABSIRACT 

An investigation is made of the wave forecasting method of Pierson, 
Neumann and James of New York University, published as lle Os Pube Noe 603. 
The ocean "swell" spectral energy.as predicted by, the application of this 
forecast method is compared to the observed spectral enorgy present as 
estimated from ocean wave recording devices installed at the surfacé and 
et a depth of 80 feet, in the coastal waters off Davenport, California. 

The comparison of the observed and forecast spectral energy for nine 
verification times taken at twelve hour intervals from 28 January 1954 to 
1 February 1954 suggested the following conclusions: 

1. For the conditions of the investigation, it was the opinion of the 
authors that the "swell" forecast was operationally useful. 5 

Ze The data of this invostigation suggested the possibility that the 
forecast energy spectrum is too high in the low and middle frequencies and 
too low in the high frequency bands for high wind speedse 

Se The data of this investigation indicated that at certain verifica= 
tion times, energy reductions due to crossesea interference were important. 

The authors wish to expross their appreciation for the advice and 
encouragement of Professor Je B. Wickham undor. whose direction this study 
was carried out. 

Appreciation is also expressed to Professor Re Le Wiegel, of the 
Institute of Engineering Research, Berkeley, California, who mado 


available the wave records upon whioh this study is based, 
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CHAPTER 1 
| INTRODUCTION 

Ihe value of reliable techniques for forecasting ocean wave proper- 
ties for seagoing and coastal operations, particularly military operations, 
is obviouse 

This problem has engaged the interest and attention of some of the 
foremost scientists of various ficlds over the years. Somo of tho contri- 
butions have been notable. 

The theorics of classical hydrodynamics were applied by Lamb ( 1] 
in 1879, which yielded important properties of oortain relatively simple 
sine-wave forms. 

In 1955, Cornish[ 2], published the results of a large number of 
carefully made visual observations of wind-waves, both in deep water and 
near the shore, which wore the first of that kind and which undoubtedly 
stimulated further work of this nature. 

The developmont of continuous wave recording devices [5] expanded ihe 
wave data available and made possible the many statistical studios which 
have been carried out. 

Sverdrup and liunk (4) published a wave forecasting method in 1947, 
for predicting the generation, decay, and transformation of the 
"significant" wind-wavos. 

Bretschneider, in 1951, published a modification forecast [5] of 
the Sverdrup and Munk technique based on more observational data. 

In 1953, Pierson, Neumann and James of New York University pub- 
lished a new theory and corrosponding forocast technique (6 | e This 
theory and method are in many respects unique and depend largely on same 


assumptions which have had limited testing. 
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Sinoe this last theory and technique has been adopted for use by the 
United States Naval Service, a study and evaluation of the mothod is the 
purpose of this papere A study of this nature presents many difficulties 
and requires data and facilities available at very few institutions. These 
facts explain the paucity of published results of such studies. To the 
knowledge of the authors, no such studies and evaluations havo been made 
previously of the Pierson, Neumann and James wave forocast mothod, hence- 
forth roferred to as the N. Y. Us method, for Pacific Ocean swell, although 
the waves forecast in a generating arca have been studied 16]. 

The study divides itself naturally into four distinct phasos: 

le A wave forecast for a given period of time at a specified location, 
using the N. Y. U. method, 

2. The analysis of wave records from a wave recording device for the 
same location and same time to obtain the actual wave conditions. 

Je The comparison of the forecast and observed wave properticse 

4. The results or conclusions to be drawn from the comparison. 

The data and facilities necessary to carry out this study are listed 
below with their sources: 

І. A forecast location at which suitable wave records have been 
kept for a convenient period of time. The site chosen was Davenport, 
California, where the University of California, Berkeloy, had a series of 
recording devices (7] installed from the period of Nov. 1952 to Dec. 1954, 
Wave records from Jan. 1954 to Feb. 1954 were made available to the 
authors for use during this study. 

де A rapid computing device for the analysis of the wave records [8 ] 
to obtain the distribution of wave enerzy as a function of frequency (9). 
The U. S. Naval Postgraduate Sohool is equipped with а 080 102-А medium- 
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speed seneral-purpose electronic computer which was available for this purpose. 
Se Synoptic surface weathcr maps for the period of and proceding 
the forecast time for the purpose of proparing the detailed wave forecast. 
Such maps were made available by the Department of Mechanical Engineering 
Wave Research Laboratory, University of California, Berkeley. 
4. A submarine contour ohart for the area adjacent to the forecast 
location. U. Se Coast and Geodetic. Original Survoy Sheet Noe 5266 of 
20 April 1935, was mado availablo by the Hopkins Marine Station, Pacific 


Grove, California. 





CHAPTER 2 
TIE THEORY AND THE METHOD 

The objective of this chapter is to discuss the portions of the 
Ne Ye Ue Method [6] and theory which relate diroctly to this invostiga- 
tions. 
2A The Spectrum of Wind=Generated Waves 

The primary basis for this forecast method is the assumed distribu- 
tion of ocean wave energy for a given wind speed as a function of wave 
frequency in a windewave generation area, or the energy spectrum This 


distribution is represonted mathematically by: 
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A 
[А = Eus ге” атар Bye (1) 


Where c is a oonstant of proportionality and the other terns are defined 
on page VI. Figure 1 shows a plot of this function for wind specds of 20, 
50, and 40 knots, and tabulated values of face}? as a function of wind 
speed and frequency are given in Appondix I. 

The form of this function was assumed by Neumann [10] based on 
theoretical considerations and empirical evidence. То а large measure, 
the acouracy of this method will be determined by how closely tho assumed 
spectrum approximates the true spectrum. Such determination will be 
possible only as increased test results become available. 

241 Energy in the Spectrum 

The integral of Equation (1) for a given wind velocity with respect 
to frequency is proportional to the energy of the included spectral com- 
ponents. This measure of spectral energy, between the frequencies Т 


and $2 , із defined as E and is given by 





2 | 
Е = ft (^ аг (2) 
Г. 


and has the dimensions, length squared. 

Well-tested statistical relationships [11 ] between E and wave heights 
have been demonstrated. These relationships are given in Table 1 and are 
used as a basis for all wave height dotorminations, E bein; tho forecasted 
parameter e 
242 Properties of the Spectrum 

The wind-gencrated sca for a given wind speed is said to be fully 
developed when all possible wave components in tho spectrum between $ = O 
and f = 02 are present with their maximum amounts of spectral energy. 

This state is reached only after the wind has blown for a sufficient 
length of time ("duration") and over a sufficiently long distance ("fetch"). 

If the wind duration or fetch is insufficient for the fully developed 
state, a lower energy level will result. The spectral components form in 
the high frequency end of the spectrum initially and the lower frequency 
components are added progressively and grow as the wind duration increases 
until full development occurs, providing the fetch is not limiting and the 
wind continues. One important feature of this theory shows the partially 
developed spectra for limiting conditions of fetch and duration for 
various wind speeds. 
2А5 The Co-cumulative Spectrum 

The energy spectrum is converted to a more convenicnt form by 
integrating it over the entire frequency range. This procedure is illus- 


trated in Figure 2. It will be noted that the ordinate of the cumulative 





integral can DS read directly in values of E, which is proportional to all 
of the wavo cnergy in spoctral components equal to or groater than the 
abscissa, fe The functional relationships betwoen energy, wind, duration 
and fetch may also be included on the cumulative spectrum. 

Figure Ó and Figuro 4 illustrate the N. Y, U. co-cunulativo spestra 
given as functions of duration snd fetch respoctivoly for verious wind 
speeds. These are the curves normally used in the N. Y., Ue methode 
2B The Wave Forecast in tho Generating Area 

The foregoing are required for forecasting the wave conditions in a 
windewave generating areae Knowing the wind spoed, the wind duration, and 
the fetch length, the E-value and the range of wave periods with import- 
ant energy may be obtained diroctly from the co-cumulative spectrum. The 
sca in the generating aroa has an irregular, short-crosted, appearance as 
the spectral components randomly amplify and damp each other. The 
spectral components also move in a varioty of diroctions, as the wind 
action is nov uniform, but of varying strength and Вара бабах 
2C The "Swell" Forecast 

When wind-generated waves loave the generating arca, a transforma- 
tion takes place in the appearance of the sea surface as tho wind no 
longer acts on tho wave form. As the various spoctral component wave 
trains leave the wind areca, the apparent waves become regular and long 
crested, and arc termed "swell." When "swell" arrives at points distant 
from the genorating area, the waves have only a fraction of the onergy 
originally present in the fetch. 

201 Assumed Processes Involved in Energy Reductions 
Pierson, Neumann and James account for energy reductions at points 


distant from the fetch by two processes, namely 


6 





integral can be read diroctly in values of E, which is proportional to all 
of the wave energy in spectral components equal to or groater than the 
abscissa, f. The functional relationships between energy, wind, duration 
and fetch may also be included on the cumulative spectrum. 

Figure 3 and Figure 4 illustrate the П. Y. U. co-cumulativo spectra 
given as functions of duration and fetch respoctively for various wind 
speeds. These aro the curves normally used in the HN. Ye Ue mothod. 
2B The Wave Forecast in the Generating Area 

The foregoing are required for forecasting tho wave conditions in a 
windewavo generating area, Knowing the wind spoed, the wind duration, and 
the fetch length, tho E-value and the range of wave periods with import- 
ant energy may be obtained directly fram the co-cunulativo spectrum. The 
sea in the generating area has an irregular, short-crested, appearance as 
the spectral components randomly amplify and damp each other. The 
spectral components also move in a variety of directions, as the wind 
action is not uniform, but of varying strength and 22 он, 

22 The "Swell" Forecast 

When wind-generated waves leave the generating arca, a transforma- 
tion takes place in the appearance of the sea surface as the wind no 
longer acts on the wave form. As the various spectral component wave 
trains leave the wind arca, the apparont waves becomo regular and long 
crested, and are termed "swell." When "swell" arrives at points distant 
from the generating area, the waves have only a fraction of the onergy 
originally present in the fetch. 

251 Assumed Processes Involved in Energy Reductions 
Picrson, Neumann and James account for onergy reductions at points 


distant from the fetch by two processos, namely 
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l. Dispersion 

2. Angular Sproading 
252 Dispersion 

When the spectral components leave a generating area, the associated 
energy front of each travels at different group’ velocity. The group 
velocity for any spectral frequency is given by 

CT Т | | (5) 
Dispersion, then, is defined to be "the spreading-out effect caused by the 
different group velocitios of the spectral frequencies in the original 
disturbance at the source." [6] Therefore, to determine for a given time 
the total spectral energy at a specified distance from the fetch, based on 
the effects of dispersion, it is necessary to compute the range of fre- 
quencies which could be present at that time. This may be accomplished by 
considering the distance the enorgy must travel, the group velocity, and 
the time the particular spectral component was formed as a function of 
wind speed, duration, and fotch. 

The N. Y. 0, method includes some basic "filters" which are ossen- 
tially pre-arranged formulations of these variables, for cortain recurring 
situations. The use of these "filters" permits the calculation of a Tr, 
upper period, and Tr, lower periode Ty and T; dofine theo upper and lower 
limits of the spectral energy present at the forecast time and location, 
for & given wind speed, 

203 Angular Spreading 

Ihe considerations of paragraph 202 do not take into account the 
effect of variability in the direction of wave motion and the fact that 
the waves are short-crested. The component wave trains vary not only in 
frequency but also in direction within tho area of generation. Thus, as 
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the spectral components leave the generating area, they spread out 
angularly as well as "disperse," 

The assumod distribution of relative wave energy as a function of 
direction is illustrated in the lowcr curve of Figure 5. Since the area 
under this curve is unity, when integrated cumulatively over direction it 
gives that fraction of the energy to be found between any two directions 
(measured from that of the mean wind). The ordinate of the upper curve in 
Figure 5 is then expressable in percent. Dispersion and angular spreading 
computations are illustrated in chapter 3.7 

Figure 6 illustrates the measurement of the angles 9, and Өз, which 
are the parameters used in conjunction with the upper curve of Figure 5 
for evaluating the effect of angular spreading. Fhe convention adopted 
for measuring @, and 9; is that clookwise angles are positive and counter- 
clockwise angles are negativee | 1 

The difference between the percentages obtained by entering the upper 
curve of Figure 5 with 6, and 8z is the percentage within that range of 
angles of the spectral energy computed from the dispersion envelope and 
the energy spectrum. 

Тһе bell-shaped distribution of energy with direction, ligure 5, was 
inferred from limited evidence [12] and to the kmowledge of the authors, 
its exact form has had limited observational verificatione 
204 Other Sources of Possible Energy Loss 

Pierson, Neumann and James [6] have stated that in many cases, 
dispersion and angular spreading are the only mechanisms necessary to 
account for the energy decrease observed in wave spectra at locations 


. remote from the generating area. They held that there may be important 





energy lost due to the effoct of cross-sea interference," but havo 
expressed the opinion that losses duc to viscosity are generally negli- 
sible. The effects of cross and opposing winds on wind-penorated swell 
arc not known, but undoubtedly would modify the spectrum of the swell to 
an extont dependent on the wind speed and the spectral composition of the 
cross-sea. It is believed that at least a qualitative evaluation of these 
effects will be possible as the volume of test results become greater. 

Ihe direct measurement of attenuation presents serious difficulties. 
However, a careful analysis of the actual spectral energy at a wave 
recorder as compared to the forecast spectral energy may indicate the 
presence of processes which cause enersy losses. For instance, theory 
indicates that viscosity effects, if they are important at all, reduce the 
energy primarily in the high frequency end of the speotrum [8]. If high 
frequency energy from a distant fetch is forecast to errive at a wave 
rocorder at a given time, and fails to arrive, it might be inferred that 
viscosity effects were responsible for the energy loss. This possibility 
is discussed further in Chapter 6. 

205 Goals of the Investigation 

The goal of this paper is to make available the results of the 

comparison of the spectral energy forecast by the N. Ye Ue mothod and the 


spectral cnergy as obtained from a wave-recording devioe, 


"The assumed effect of cross-sea interference is that tho swell with 
spectral components of the same frequency as the cross-sca wave trains 
will form unstable interference patterns with the latter. This will re- 
sult in attenuation of the swell having speotral energy in those frequencies. 





It was intended that these results would permit conolusions parti- 
cularly in the following areas: 

lo The assumed form of the energy Spectrum, 

2» The assumed form of the wavc-energy directional spectrum. 

e The effects of viscosity and oross=seas on the speotrum of the 
swoll waves. 

4. A general vorification of the system which implicitly includes 


the items above. 
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CHAPTER 3 
THE FORECAST OF THE ENERGY SPECTRUM 

This chapter contains a brief summary of the procedures followed by 
the authors in proparing a forecast of the energy spectrum for a given 
observation time, 

When a breeze comes up, the sca surface instantaneously becomos 
covered with tiny ripples which form more or less regular arcs of long 
radius, They increase rapidly in height as energy is transferrod to the 
waves by the pushing and dragging forces of the wind. бооп a stato of 
irregular composite wave motion, called the sea, exists. The fully 
developed state over a long fetch with a sufficiently long duration 
finally will be attained when the wave-generating tractions have increased 
tho total wave energy to such a point that dissipation balances the work 
done by pushing and dragging forces on the sea surface. The — amount 
of wave energy accumulated in the composite wave motion is now distributed 
over a wide range of frequencies. The curvos and assumptions of the 
Ne У. Ч. method were used in forecasting the energy spectra in the fetch 
areas He Oe 605 contains several filters that can be applicd to account 
for the reduction in the energy spectrum as its components move from the 
fetch area to an observation point. The authors, after considerable 
investigation and study, oonoluded that the idealized filters in H. 0. 


605 would not be oonvenient for use in this study. To achieve the de- 


sired results, a new method of solving for dispersion and angular spread 
ing was devised. The new filter used when solving for the effects of 
disporsion allowed for wind variations in the fetch and resulted in a 
continuous dispersion envelope of periods present at the forccast point. 
The angular spreading factor allowed for any orientation of the fetch. 


Ti 





ЗА Selection of Data 

Appendix II describes the elaborate instrumentation used by the 
University of California to record continuously the water level at 
Davenport, California. An examination of the wave records available 
showed that the average wave heights were higher in the winter at Daven- 
port due to well-developed storms in the North Pacific area. Therefore, 
the months of January, February, and March 1954 wore selected for closer 
studye 

Weather maps at six-hour intervals were then studied for these three 
months, All possible wind circulations that resulted in fetches acting 
with components toward Davenport were recorded. An envelopo showing the 
periods (Ty and T;,) and their estimated times of arrival (ETA) was com- 
puted for each fetch. These data (Ty, Ty, ETA) were tabled and plotted. 
After detailed examination, the period between 06502, 25 January 1954 and 
04502, 1 February 1954 was selected for thorough investigation and study. 

This time interval was the most favorable for the following reasons: 

le A stationary nearly uniform fetch prevailed. It is discussed in 
Appendix III and illustrated in Figures 7 and 8. 

2e Representative wave records were available in "fast timo"? with. 
a time interval of twelve hours, as explained in Appendix II. 

5. Fetches A and B? contained weather station Blip. Y Po ls 
manned by experienced aerographers. Therefore, its woather reports were 
assumed acourate. "Р" also reports sea-water temperature (Tg) which 


lrigure 9A is en example of this fram 06302, 25 January 1954 to 04502, 
l February 1954. 


é"Fast Time" is used to indicate the portions of the record where sea 
level was recorded on the scale of 6 inches per 1 minute. 


3 See Appendix III. 
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could be used in computing the stability factor, which is used in 
estimating the wind speed. 

4. Only fetohes A, B, and C, as defined in Appendix III and illus- 
trated in Figures 7 and 8, were believed to be producing significant 
energy at Davenport e 

5, Looal winds in the vicinity of Davenport were blowing from the 
land toward the sea and were light, 

6, All instruments were cleaned and calibrated on 24 January 1954 
and assumed to bo in good working order, 

Considerations 1, 2, and 4 established wave forecast times at 
Davenport which are as follows:. 

$] 04502 28 Jan 1954 
$2 16302 28 Jan 1954 
tz 04502 29 Jan 1954 
t, 16502 29 Jan 1954 
tg 04502 50 Jan 1954 
tg 16502 30 Jan 1954 
ty 04302 51 Jan 1954 
$3 16502 51 Jan 1954 
бо 04502 1 Feb 1954 

Hereafter, the forecast time will be referred to as t;e 
óB Compilation of раба 

Table 2 reoords the data measured from the weather charts illustrated 


in Figure 7. 


lsee Table b, 
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Table 3 records the data computed from parameters listed in Table 2 
by the following proceduros (symbols are all defined in Tables 2 and 3): 

1. Stability Factor (V,,/V,) depends on Tg minus T, and is determined 
from Table 5. 

2. Isobar curvature correction (1,,) depends on isobar curvature (I, ) 
and stability factor and is determined from Table 6, 

3. Geostrophic wind reduced to sea surface D 

Veg = (Vg) (Igo) (stability factor) (4) 

de Surface wind speed (V) obtained by comparing Уор, the reported 
surface wind, with Voge 

5. Limiting fetch (F, ), upper period C and E were read directly 
from Figure 4 using Ve 

6. Equivalent duration (Tp) was established by camparing generation 


time of last V to present V, for an equal E-value. 


Te Travel time 225) 


Т - Ко Ра 7 37 | (5) 
ob = 
21 515 TU 


Where Ty leaves the wind area from near the windward side of the fetch. 


~ 


8. Lower period (Ty) 


R 
тг = er. (6) 
1.515 Ton 


This T, will have the same ETA as Ty because it leit from the leeward side 
of the fetch. 
9. Duration time (T,) 
Та Е ТЕ HAT (7) 


Where AT 3 6 hours, or map interval. 


“Бөө Figure 7 
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10. Estimated time of arrival (ETA) 
ETA = map time < Tob (8) 

90 Development of the Filter 

The values of El must be reduced by filters as the spectral components 
move from the generating area to the obsorvation point. Chapter 2 explained 
the theory of filters involved in the N. Y, Us method. Two processes act 
to limit the energy at the observation point: dispersion and angular sproad- 
inge 
D Dispersion 

Dispersion results from the presence of different group velocities 
among the spectral components in the fetch. The writers, to accurately 
account for the effects of dispersion, plotted Ty and T, аза function of 
time.^ Ty left from the windward side of the fetch with R, Æ Fie -~ Fh 
| distance to travel. Ty, moved out of the leeward side with distance Rọ to 
travel. Only periods between T, and T, were theoretically present at 
Davenport because of the dispersion influence, 

Fetches A, B, and C were fully=developed wave generating arcas; 
consoquently, upper period waves (Ту) were also leaving from the leeward 


side with travel distance Ryo Their ETA was computed by 


R 


Top s 9 
pueris T. P 


and Equation (8). Due to a shorter travel distance this T, arrived ahead 


of a similar Ty which left from fetch windward. 


looo Table 3 
“Soe Table 9A 
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To include the first arriving waves from the generating area, the 
waves created as the wind continued and the waves leaving the fetch area 
when the wind ES Figure 9A was modified into Figure 98, а dispersion 
envelope. А dispersion envelope is a plot of Ty and Tz, present from the 
same fetch as a function of the same ETA. The DIN ODORE plotted from 
values obtained through the use of Figure 4 and Equations (5), (6), (8), 
and (9). The results are shown in Figure 9B which illustrates the various 
limiting envelopes of periods due to dispersion affecting Davenport before 
and during the forecast periods from fetches А, В, and Ce 

The resulting reduction in E due to dispersion is illustrated in 
Figure 10 for t) and Figure 11 for ty. 
óE Angular Spreading 

Angular spreading 2 from the variability in wave direction in 
апа after leaving the fetch area. The angular spreading factor is related 
to the orientation of the observation point with respect to the generating 
area and is therefore dependent upon the parameters 957, ӛлі; Әзу, and 
m which were measured from the synoptic weather chart. lor times near 
the end of a fetch life, important wave energy may come from the rear of 
the fetch (fetch windward). At these times the angular spreading factor 
decreases in value, consequently, less encrgy will be forecast to arrive. 

The factors for the fotch leeward and fetch windward were obtained 
from an angular spreading factor graph." The difference between these 
two factors was (from leeward to windward) as great as 203. The authors 


considered the possible error introduced by using an angular spreading 


lsee Figure 7 
“See Figure 5 
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factor averaged over the entire fotch as unacceptable. To accurately 
determine the angular spreading value as a function of distance from the 
fotch front, Equation (10) was used. 
А 1,515 (Top) T = Ro (10) 

“AD is the distance in nautical miles measured from the leeward side of 
the fetch into the generating m 
Top is the travol time. 
T is any period between Ty and T; for t; as determined by the dispersion 
envelope e 
Rg is the distance from observation point to the center of fetch leewarde 

In Equation (10), T was replaced by l/frequency, a 4D was computed 
for each 0.005 chango in frequency for each ti involving fetches A, B, and 
Ce The distance, AD was thon caloulated and plotted, and a corresponding 
9- and 6, measured. These angles, applied to the angular spreading graph, 
yielded the angular spreading factor. This procedure was used when it was 
evident that the waves were not leaving from the front of the fotche 

The resulting reduction in E due to angular spreading is shown in 
Figure 10 for tj and Figure 11 for ty. 

It is considered significant that the use of Equation (10) allowed 
for variations in angular spreading due to variations in (Ro FAD). 
Figure 12 is a graph of angular spreading factor with respect toAD, for 
fetch A. This factor did not decrease at a constant rate in going from 
fetch leeward to fetch windward. Graphs, not included hero, yielded 
similar results for fetches 5 and C. These variations in rate of change 
of the angular spreading factor depend upon the oricntation of the various 


fetches with respect to Davenport. 


lsee Figure 7 
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SF Shallow Water Effects 

Shallow water effects of refraction, shoaling, percolation, dif- 
fraction, and friction influence E as the deep water wave energy approaches 
the Davenport area. ` 

Refraction was determined using the method presented in IH. 0. 254, 
end a wave refraction diagram constructed using 0, 5, Coast and Geodetic 
Original Survey Sheet Noe 5266 of 25 April 1935. The chart shows tho 
bottom contours seaward of Davenport to be nearly straight and parallel. 
The range of deep water wave directions was very small, therefore refraction 
was treated as if all deep water waves arrived from one angle. Shoaling 
effects were also obtained from Н. 0. 234, Plate I. 

Percolation and friction were assumed negligible due to the narrow 
width of the continental shelf. Diffraction was assumed to be negligible. 
All instruments were locatod in open water well away from rock formations 
or man-made structures e 

The shallow-water effeots of refraction and shoaling wore combined 
into one correction factores Their influence on E is shown in Figure 10 
for tj and Figure ll for tq. Figures 10 and 11 show the resultant ocean | 
wave spectral energy forecast to arrive at the recording instruments for 
9) and t4e The area under the curve can be integrated by a planimeter 
to evaluate the total energy present for any forecast time. The results 
are shorm in ohapter 5, and are plotted with respeot to time in 


Figuro 14, 
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CHAPTER 4 
SPECTRUM ANALYSIS 

Ihe objective of this chapter is to discuss tho general procedure 
and theory employed in obtaining the energy spectrum from the record of a 
wavo=recording devices ў 

In this study, wave records were used from both a surface wave гесога- 
ing device and a pressure recording device." The distinction between the 
two is that the former records surface wavos and the latter records wave 
pressures at the depth in which it is installed, 
4A Theory 

Many of the features of surface, and below surface waves have been 
ined by Pierson and Marks 18) and Neumann [15] basod on the assumption 
that a wave record as a function of time is the result of a Gaussian 
process. The theory was developed by Tukey and Hamming [9] among others, 
In the publication cited, the latter oxplain a method which may be used 
for estimating the power spectrum of a steady-state sea surface using the 
Fourier transform of the autocorrelation function. A text which treats 
correlation functions and power spectra in detail is Information Theory[i7?} 
by Goldman. 

Timme and Stinson n4] present ourves of the pressure response 
factor as a function of wave period, T, and depth, He The pressure 
response ‚factor is a moasure of the attenuation of wave enorgy with depth, 
4А] iumerical Estimation of Non-normalized Autocorrelation Function 

and the Pressure Power Spectrum 


The results of Tukey end Hamming Г9) are used to derive the power 


& 
\ 


l$ee Appendix III 
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spectrum from a pressure record placed at the sea floor in a depth of 
water, lle The possible errors which enter into the numcrical analysis 
have been considered. The corrections which must be appliod to the raw 
results and the degree of reliability have boen determined, 

The formulas for numcrical estimation of the powcr spectra as 


presented by Tukey and Hamming[9] end as used by Pierson and Marks (8] are 





given by: 
to = ty, = $5 = $22 ---- Бро бі а At (11) 
цер | 
г: P(ty) P(ty.p) (P=0,1,2, ---- m) (12) 
nz 
пе] p 
L = = (Q of2 > Qp cos SPA. ¥ qm CosTh) (13) 
p=] 
й (h=0,1,---n) 
Un = „2510 2054 Ly 4023 yyy (hrO,1, ---m) (14) 
Ау з г. (15) 
6° =. FE (16) 


Where = хары 

Equation (11) states that li-1 equally spaced time intorvals aro 
marked off on the record by points spaced At units apart. The values 
of pressure at these points in terms of departures from the mean are then 
tabulated as P(tj), P(t2) ---- P(t,), ---- P(t);), N numbers resulting. 

The non-normalizod lag coefficients are then found by the use of 
Equation (12). Qp ranges from Qo to Qu» or m 7 1 total numborse Qo is 
the Eph max, or the total spectral energy uncorrected for the scale factor 
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of the record, at the depth H. 
The sumnation of Equation (12) was acoomplished by the use of a 


CRC-102-A computer.! 
ТГ(в-2) , 
b c 
AA лай are found fran Equation (13). A plot of the "raw" estimates 


e "raw" estimates of power in theJ4 band frmu = 


with values of h result in a quito jagged graph. The jaggedness is intro- 
duced by the measurement of P at disorete points and by the finite length 
of the record. A smoothing device is given by Equation (14). 

Ihe final number, Uu, is an estimate of the mean value of the power 
system from 1 == фо TE. At the point given by Equation (15) 
the best estimate of the value of EN (797); is givon by Equation (1252 
The result is 8 smoothing of the "jazged" curvo of Equation (15), and this 
continuous curve is the best estimate of the pressure power spectrum. 

For the purposes of this investigation, the computcr program for the 
Fourier analysis was limited to Equation (15). The refinements, Equations 
(14), (15), (16) were not included in the solution, therefore, the resulting 
Spectra are the "raw" or unsmoothed curves. This limitation was dictated 
by time considerations as the preparation of such a program is extremely 
longthy. These refinements will be included in the computer program as 
time permits and will be available for similar investigations in the 
future, 
` 4A2 Reliability of Results 

The final results of paragraph 4Al are estimates. Estimates can be 
in error, and the important part of this analytical method is that it 


also presents an ostimate of the error. Tukey and Haming [9] have shovm 


lone computer programming of the equations used in this investigation 
was accomplished by Professors Je Be Wickham and Н. М. Martinez of the 
0, 5, Naval Postgraduate School. 
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! 
that tho errors follow a Z? distribution with f degroes of freedom given 


by 


N ->m 
=. Sern е. 


fs 2 (17) 
Givon f, thon the use of Table (7) yiclds information on tho reliability 
of the estimate for each frequency band obtained by the analysis. 

For 20 degrees of freedom, theU, obtained from Equation (16) may be 
multiplied by 1.8 and by.65. The true value of U in the band under analysis 
will lio between 1.8Uy and „65у n 905 of the time. 
4B Procedure 

For this investigation, the verification periods were of the order of 
twenty minutes on the wave rocord. The +t was chosen as two seconds as 
determined by the criterion of Pierson and Marks [8]. The total number 
of measurements per wave record, N, was of the order of 600, or pro- 
gremming convenience, the total number of lags, m, was taken as 43.1 

The total number of consecutive verifications, 12 hours apart, was 
nine.” It was originally intended that an autocorrelation program and a 
Fourier analysis would be carricd out for each of these times, However, 
somo initial difficultics in the programming routine rostricted the time 
available for the investigation, and time did not permit the fulfillment 


of this objective, 


Im future work, the wave record should be analyzed over a longer time 
interval resulting in greater value of ll. This will increase the degrees 
or freedom and permit better resolution of the spectral curvee The 
number of lags taken should be increased, 


“Seo pp 15, Chapter 3. 
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CHAPTER 5 
STATELENT OF RESULTS 

The purpose оҒ this chapter is to present the results of this 
investigation. 
5А The "Filter" 

It is the opinion of the authors that the "filter" devised for the 
energy spectra forecast in this study represents an accurate, practical 
method of solution. This is particularly true for the case of marked 
variation of wind speed with time in a generating areae Гог example, in 
the case of increasing wind speed with time, energy fronts generated by 
tho earlier, lower speed winds will frequently be overtaken by energy 
fronts gencrated by the later, higher speed winds. This additional 
energy, from the higher speed wind might be overlooked if a plot of the 
"dispersion envelope" with time were not maintained.” 

The filters and examples desoribed in the N. Y. Ue method 6 
primarily treat the simplified cases of constant or slowly changing wind 
speed. This permits the energy to be obtained directly from the 0.0.5. 
curves for one wind spcede 

It must be understood that the "period envelope" must be identified 
precisely with the wind speed which generated it. It appears that the 
filter herein described? adequately handles the moro complicated 
situationse 

It appears that in forecasting by this method, using the filter 
presented herein, the forecaster will frequently need to work directly 


with the spectral curves in addition to the C.U.S. curvese For this 


!Seo Figure 9B 
“See Chapter 3 
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reason, the tabled solution to Equation (1) is inoluded in Appendix I. 
It is recommended that personnel using this method plot the curves of 
the spectra for various wind speeds on an enlarged scale for tho purpose 
of representing the spectrum present at the locality and time of interest. 
5B Verification of the Forecast 

Spectrum forecasts were prepared for Davenport, California, for the 
times ty through ty as described in Chapter 5. These forecasts were 
prepared for two geographical sites? having water depths of 50 and 81.5 
feet respectively. A description of these sites and the instrumentation 
of each is presented in Appendix II, 

Data available for verification of these forecasts were the 
following: 

1, Significant wave height measurements for the sea surface,water 
depth 50 feet, hereafter called site I. These data were computed from 
sea surface wave records by personnel of the University of California, 
Berkeley, and made available to the authors for this study. Significant 
wave height, H 1, is defined as tho average of the heights of the one- 
third К ОООО оз. E values may be reoovered by the statistical 
relationship, 

й = = 2,88 VE | (18) 
A plot of the E values obtained in this manner for the times t] through 
to is illustrated in Figure 14. 

Ze A sea-level wave record for site I obtained from a Beach Erosion 

Board, Step-resistor gauge as described in Appendix II. These records 


were analyzed by the procedures described in Chapter 4, and an estimate 


15оө Figure 13. 
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of the E value for timos tj and to were obtained, 

3. A "pressure" record for the site in the water depth of 81.5 feet. 
Ihis site will be referred to as site II. This pressure rocord was ob- 
tained from a Mark IX pressure-head instrument described in Appendix II, 
The instrument was installed 1% feet from the bottom at site II or in a 
water depth of approximately 80 feet. These records пого analyzed by the 
procedures described in chapter 4 and an estimate of the E-value for 
times tz through tg were obtained. 

As described above, the "reoords" available were obtained from two 
different instruments. For the times t} and tə, the surface step- 
resistance device was available at site I. For the times tz through $9 
the pressure=head device, installed at a depth of 80 feet, was available 
at site Il. 

Figure 14 is a plot of E-valuos for the times ty through tg for the 
sea surface at site I, as estimated by: 

le The forecast. 

2. The H 1 computations provided by the University of California, 
as described Eo The 90% confidence limits [6] are plotted at each 
pointe 

96 The writers! analysis of the water level record of the step» 
resistance gauge as described above. The 90% confidence limits for 
these points are also included. 

Figure 15 is a plot of E-values for the times t; through tg for 
site II for a depth of 80 feet, as estimated by: 


le The forecast. 


15ее рр 22, Chaptor 4. 
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2. The writers'! analysis of the "pressure" record of the instrument 
as described above. The 907 confidence limits for these points are also 
included. 

The procedure for obtaining the forecast spectra and b-values for 
verification times tz through tg at site II was the following. 

A surface spectrum forocast for a given verification was prepared at site 
11, This spectrum was reduced to a depth of 80 feet using the pressure 
response attenuation factors of Timme and Stinson ÜŞ] as described in 
Chapter 4. An example of the results of such a procedure is illustrated 
in Figures 16 and 17 for verification times ts and tg respectively. 

As stated earlier, the original intention was to obtain the computed 
spectrum for each verification time. However, because of some initial 
diffioulties oncountered in the computer programs disoussod in Chapter 4, 
the time allotted for the use of the computer did not pormit the fulfill- 
ment of this objective. 

A disoussion of the figures described above, and the conclusions of 


the authors will be presented in Chapter 6. 


15се pp 22, Chapter 4. 
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CHAPTER 6 
d CONCLUSIONS 

The purpose of this chapter is to compare tho obsorvations DT the 
forecasts and to suggest possible explanations for any differenoes noted, 
64 General Comments on the Results 

An examination of Figures 14 and 15, which compare forecast and ob- 
served energy, and Figure 22, which compares the forocast and observed 
significant height, indioates certain interesting featurose 

1, It is noted that the general agreement is quite goode The fore= 
cast energy peaks occurred at the time predicted and the energy was of the 
correct order of magnitudo, 

де The forecast energy was in general equal to or greater than the 
observed, 

Se At the point of greatest difference the foreoast significant 
height was about one and one-half feet greater than the observed. This 
would indicate that for these conditions, the forecast would be of 
Operational value. 

4» It is noted that the points of closest agreement occur where the 
observed energy is at or near a peak. 


6B Possible Explanations for tho Differences between the Observed and 
Forecast Energy . 


For the purpose of suggesting possible explanations of diffcrences 
between the observed and forecast energy, the following items are con= 
sidered to be possible sources of error: 

le The theory. The foreoast energy depends upon 

ае the energy spectrum and, 


be the directional spectrum. 
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If the assumed form of eithor of these curves is incorrect, tne onergy 
forecast would probably be in error. 

2e The application. In applying the methed, it is unavoidable that 
the selection of data is largely subjective. It is possible therefore that 
sizable errors could bo introduced. Errors in the following might be 
particularly critical: 

ае geostrophic wind analysis 
be surface wind calculation 
Ge measurement of fotch parameters. 

Зе Energy transformation processes. Processos or conditions which 
are believed to affect ocean swell wave energy and are not accounted for in 
the forecast are the following: 

ае cross=sca interference 
be cross or opposing wind 
Се viscosity effects. 

4. Refraction and shoaling effects. 
6C The Energy Porocast 

Figure 14 compares the forecast value of E with the observed at site I 
at the sea surface for the verification times ty through tge It is noted 
that the forecast E value lies above the 90% confidence limits of the 
‚ observed E value at the times Be tas ба» and tge The inference might be 
dram that at these timos onergy loss processes, not considered in the fore- 


cast, are of importance. It is believed significant that a well defined 


"Additional computations of these data indicate that the correct value of 
E is approximately identical with the forecast value. 
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aroa of crossesea was present at such location and timos as to have inter- 
ferod with the "swell" arriving at the times tz, t,, and tg from fetches A 
and В (Appendix III). Crosseseas aro believed to cause the summation of the 
spectral components to become non-linear giving rise to the convorsion of 
wave energy to turbulent onergy by formation of whitoecapse Figure 18 
illustrates the looation ond orientation of the cross=sca relative to the 
generating fetches and the forecast sites at Davenport. Figuro 19 displays 
the spectrum forecast for the cross=soa generating area at six hour inter- 
vals over the period when it is most clearly defined. For reasons vi 
stated, computed spectra for the verification times under discussion, аге 
not available. It is believed that suoh spectra would havo aided greatly in 
the determination of the portion of the spectrum affeoted by the cross-sca 
interference. However, it was noted that if the foreoast swell cnergy were 
totally erased in the frequency band of the cross-seas, energy reductions 
of the order of those observed in figure 14 would have occurred, 

Figure 15 compares the forecast values of E with the observed values 
at site II, at a depth of 80 fect, for the times tz through бо» 

By comparing the energy difforences which ocour at tho surface 
(Figure 14) with the energy difforences which occur at depth (Pigure 15) 
it appears it might be permissible to draw inferences about the frequency 
ranges in which these differenoes occur. These inferences aro based on 
the fact that the middle and lower frequency energy would bo less 
attenuated with depth by the hydro=-dynamio filter than would the high 
frequency energy. Accordingly, if the relative error of the energy differ- 


enees at the surface was of the same order as the relative error at depth, 


it might be inferred that the energy differences were primarily in the 
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middle and low frequencies. By similar reasoning, if tho relative error 
at depth was significantly less than the relative error at the surface, 
it might bo inferred that the energy difforencos were primarily in tho high 
frequency end of the spectrum. If tho relative error at depth is greater 
than the relative error at tho sca surface the possible causes aro less 
obviouse A possible oxplanation might be that too much middle or low 
frequency energy and too little high frequency energy is forecast, the 
errors cancelling and giving close agreement with the observed energy at 
the surface. Then, as tho forecast spoctrum is reduced by the hydro- 
dynemic filtor, at depth the surplus middle and low frequenoy energy which 
was forecast would be relatively less attenuated with depth than the high 
frequency energy. Then at depth, the forecast energy by virtue of the 
surplus middle frequency energy would show a groater relative orror than 
at the surface e 

It should be notod that it is quite possible that minor vind systoms 
such as sea-breeze circulations, may have developed between the six hour 
map intervals of this study which could have developed enorgy sources not 
taken into account. Further, because site II is some distance semvard of 
site I, the total energy at site II is greater as the shoaling effect! is 
less pronounced. Since, in the above discussion, the onergy difforonces 
are compared between site I at the sea surface and site II at depth, the 
relativo error as observed from Figure 15 will ba somewhat greater than 
would be the relative error at depth at site I. 

xamining the relative errors of energy differences at the soa- 
surface and at depth for the timos tz through to, and applying the 

1 differences betweon site I and sito lI wore assumod 

negligible. 
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roasoning suggested above, the following comments are offered: 

le At the timos tz and ts the relative erros at the sea-surfaco and 
at depth are about the same. At tho time t, the relativo error is much loss 
at depth than at the surface. vollowing tho roasoning suggested above it 
might be inferred that at times tz and ts important energy differoncos 
oocurred in the middle frequencies and that at timo t, the energy difference 
was primarily in the high frequency range. As noted previously, it is 
believed that cross-sea interference acted on the "swell" arriving at these 
times. These comparisons give rise to the possible oonclusion that the 
forecast spectrum ovorprediots in tho low апа middlo frequencies and under- 
predicts in the high frequenciese 

2. At the times te and 57, the relative error of the energy diffor- 
ences was significantly greater at depth than at the surface. As suggested 
above, a possible explanation for this may bo an overforecast of energy in 
the middle and low frequencies and on underforecast of energy in tho high 
frequency band. Ib may be noted that these times are at or near an energy 
peak. It might be suggested than that for higher wind ere the assumed 
enorgy spectrum overprediots in the middle frequencies and underpredicts 
in the higher frequencies. At least such an assumption is consistont with 
these datae 

Se At the times tg and tg the relative orror at the surface is signi- 
ficantly larger algebraically than at depth. Following the line of roason- 
ing used for time ty, it might be inferred that important energy differences 
occurred in the high frequency band. The oauso of this assumed high 
frequency energy reduction is not readily apparent. Il is noted, however, 


that at both time t, and tg tho observed enorgy is at a minimum. A 
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possiblo explanation might be that at lower wind speods the assumed energy 
КОО олау overprediot in the high frequency end of the spectrum, 
60 | Sumnary 

In summary, it is the opinion of the authors that the forecast method 
for swell is operationally useful for fotches of the type active in this 
studye based on this investigation, it appears that the curves and pro- 
cedures of the method [6] , when properly applied, will result in a wave 
forecast of useful accuracye Figure 22, a plot of the foreoast significant 
wave height compared to the observed illustrates this conclusion. As noted 
in paragraph 6A above, the results are consistont with the assumption that 
onerpy reduction processes not considered in the forecast are at timos 
important. In particular, it is believed that cross-sea interforonce was 
responsible for significant energy reductions at oertain verification timos 
in this investigation. It is believed that a qualitative correction would 
be possible by reducing the forecast spectrum by a fraction of tho cross=sea 
spectrum developed. 

The data are also consistent with tho premise that, for increasing 
wind speeds, the forecast energy spectrum is too high in the low and middle 
frequencies and possibly too low in the high frequency band, and that for 
lower wind speods, an opposite trend appears possible. 
6E Future Research 

When certain modifications in the computer progrems discussed in 
Chapter 4 are completed, and accurato observational spectra are available, 
more definito conclusions will be possible, It is believed that cross-sca 
effects may become quite clear from a comparison of the east and the 


observed spectra, It is the opinion of the authors that extensions of the 
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mothods used in this investigation can provide tho quality and quantity 
of data necessary for tho ovaluation of forocast mothods and tho offects 


of аббопцаўіопе 
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APPENDIX I 
Tho purpose of this appendix is to tabulato values of [A EJE in 


Equation (1), Chapter 1, as a function of wind speed and frequency. 


2 
2 С К сс coe 
раў. (ат)бг8 ame r 


where C = 4,8 x 10% cm? sec”? 


A 
7 х 10- “4.485 х 10" 
5 P 5 ае A 
SE [A (ӘЛІ; = я 6 2 5 ү2 (2 (1а) 


The values of the [e for the solution of Equation (la) are given in 


Table (8). 
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APPENDIX II 


DESCRIPTION OF LOCALE AND FACILITIZS AT 
IHE FORECAST SITE, DAVENPORT, CALIFORNIA 


The University of California, with the pormission of Santa Cruz 
Portland Cement Company, planned, established, and maintained an ocean 
wave recording program at Davenport, California. As shown in Figuro 15, 
the pier extended twenty-two hundred feet seaward of tho beach and was 
exposed to the large waves from the North Pacific generating areas. 
Elaborate instrumentation and rocording devices were used. Continuous 
ocean wave records from 8 Nov 1952 to 1 Dec 1954 were obtainod from the 
Beach Erosion Board (BEB), step=resistor gages (parallel type) [15] and 
Mark IX ГЭ) pressure heads. The locations of these two instruments are 
shown in Figure 15. The BEB was located at the seaward end of the pier. 
Depth immediately below the BEB was 50 feet. The Mark IX was located in 
81% feet of water (below MLL). 

Fluctuations of the sea surface were transferred by electrical cir- 
ouitry to a recording center located on shore, Here the fluctuations were 
permanently recorded on Esterlino=-Angus Continuous Record Charts. Two 
time scales were used: - 

6" per 1 minute, called "fast time" 

6" per 1 hour, called "slow time" 
Fast time was recorded every twelve hours for twenty minutes beginning at 
04502 and 16302 of each daye t; (forecast time)was determined by the fast 
time records. 

A continual maintenance schedule was kept, which involved mainly 


cleaning and re-calibration. 


15ее Figure 21 
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APPENDIX III 


DESCRIPTION Or THE SYNOPTIC SITUATION 
DURING AND PRECEDING THE FORECAST T IE 


As explained in Chapter 5, the forecast time began at 0450Z, 28 Jan 
1954 and ended at 04502, 1 гер 1954. 

The synoptio situation affecting the forecast locality during and 
preceding the forecast time remained basically stationary. A low pressure 
conter (average 1000 MB) existed over the Gulf of Alaska or British 
Columbia from 12302, 22 Jan 1954 to 1250Z, 28 Jan 1954, A strong high 
pressure center (1040 MB) maintained itself over Alaska and the Bering Sea 
with a ridge extending southeastward down to 55 N latitude. Figure 7 
illustrates the general synoptic situation that existed before and remained 
nearly constant during the forecast period. 

The large cyclonic circulation in the Gulf of Alaska created a fetch 
of length (Fra) 13520 nautical miles.  Fetoh width (Ру) was 520 nautical 
miles with the distance from the center of the leeward side to Davenport 
(R » being 420 nautical miles. This fotch hereaftor will be referred to 
as Fetch "А." Fetch A existed fully devoloped from 00502, 23 Jan 1954 to 
06502, 26 Jan 1954. 

At approximately 06302, 26 Jan 1954, due to a slight reorientation of 
the cyclonic circulation, Fr. became 800 nautical miles with a correspond- 
ing increase in Ry to 960 nautical miles. This change resulted in a 
redesignation to Fetch "В." All other geographic variables та 
constant. Fetch B existed from 06302, 26 Jan 1954 to 12502, 28 Jan 1954. 

Data for fetches A and B were recorded in Table 2, Chaptor 5, and 
the fetches are illustrated in Figure 7. 


Wind is one of the most important variables in forecasting an energy 
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Spectrum. The general wind arca remained fairly constant but wind varied 
before and during the forecast period in fetches A and D. The variation 
of the average wind speed with time in fetches А and B is iilustrated in 
Figure 17. 

Figure 8 illustrates the general synoptic situation that existed from 
00502, 27 Jan 1954 to 1830Z, 27 Jan 1954. A fetch, Fetch С, was fitted 
to this circulation. Data resulting from this fetch are recorded in 
Table 2, Chapter 9. 

The wind field remained constant at 22 най during the existence of 
retch Ce 

Fetch A existed for 84 hours, B for 48 hours, and © for 18 hourse The 
subscript notation (A,,B,;,C,) was’usod to distinguish fetch sequence, with 
a nap interval of 6 hours. The fetch life started when i = 1, and ended 


when i 3 14, i = 8, and i = 4 respectively. 
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